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A search for signs of universal extra dimensions (UED) has been performed with the DO experiment, using 
events with two same-sign muons. The considered minimal UED model includes one extra dimension, and has a 
stable lightest Kaluza-Klein particle (LKP) which is a dark matter candidate. In the search, 7.3 fb~^ of DO 
data, collected in pp collisions at the Fermilab Tevatron collider at y's = 1.96 TeV, have been used. 



1 . Theory 

In the past century, many theories have been suggested that include extra spatial dimensions. An early attempt 
was made by Kaluza 1 to unify electromagnetism with general relativity by introducing a fifth dimension. 
This theory was extended by Klein [2 , who argued that the additional dimension is not observed because it is 
compactified, or "rolled up" . This branch of research was later abandoned, since as a classic theory it did not 
account for nuclear interactions, and it seemed to suggest the existence of a new particle that was not observed. 

Interest in theories with extra dimensions has increased again in the past decade, since they promise solutions 
to some outstanding problems in physics, and make predictions that can be tested in collider experiments 
such as the Tevatron or the LHC. One theory that might hold a solution to the hierarchy problem is the 
Randall-Sundrum model which assumes a five-dimensional spacetime. The standard model particles are 
confined to a four-dimensional brane, while gravity can propagate outside into the extra dimension, or bulk. The 
model of large extra dimensions by Arkani-Hamed, Dimopoulos and Dvali (ADD) [4, 5 shares the idea of gravity 
being able to propagate into an extra dimension, and aims to explain why it is so weak compared to other forces. 

In models with universal extra dimensions (UED) [5^ , not only gravity, but all forces are allowed to propagate 
into the bulk. The movement of particles in the extra dimension(s) is not seen as such, but since the particles 
have additional kinetic energy, they seem to be copies of standard model particles with a higher rest mass. In 
the following, the minimal UED model |7] with one compactified extra dimension will be assumed. This model 
has two parameters, the compactification scale and a cutoff scale A up to which the effective theory is 

valid. Here, A — 10000 GeV is assumed. The compactification introduces periodic boundary conditions, so 
that the momentum in the extra dimension can only take discrete values. This leads to the appearance of a 
series of excitations of SM particles, also called a Kaluza-Klein (KK) tower. Generally, only the first excitations 
(n — 1) can be seen at current collider experiments. The masses of the KK particles are given to first order 
(without radiative corrections) by = Mq + n^/i?^, where Mq is the SM mass. Since momentum is conserved 
separately in the extra dimension, so are these excitations, and one can introduce a new conserved quantity 
called KK parity. This multiplicative quantum number, similar to R parity in certain supersymmetric models, 
ensures that KK particles can only be produced in pairs, and furthermore that the lightest KK particle (LKP) is 
stable, making it a good candidate for cold dark matter. 

In the considered model, the excited photon (71) is the LKP. The next heavier particles are KK-leptons {£1), 
gauge bosons {Zi/Wi), quarks (SU(2) doublet Qi or singlet qi), and finally excited gluons (gi). 




Figure 1: (a) Production of pair of KK-excited quarks (QiQi) via Z boson, (b) Decay of a KK quark into a jet, two 
opposite sign leptons, and an LKP. The double lines indicate KK particles. 



The production mode of KK particles at the Tevatron is via pair production of KK gluons or quarks. These 
decay in cascades into standard model particles and LKPs. A typical production and subsequent decay is shown 
in Fig. 1] The Qi decays into an excited gauge boson {Zi or Wi) and a standard model quark. The excited 
boson decays further into a standard model lepton and a KK-lepton, which finally decays into another SM lepton 
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and an LKP (71). In the final state, there are up to four charged leptons, two jets, and missing transverse energy 
(^t) (from LKPs and possibly neutrinos). Since the emitted SM leptons may have a very low when the 
KK particles have similar masses, and thus might be below reconstruction threshold, only the presence of two 
leptons is demanded. To suppress SM background from Z — > /i/i, both leptons must have the same sign. In this 
analysis, the like-sign dimuon case will be considered, since muons have the highest reconstruction efficiencies of 
the charged leptons. 

Recently, searches have been conducted in the diphoton channel by DO [8J and ATLAS [9J for a related, but 
different model, in which the lightest KK particle 71 may decay into a graviton and a photon. However, for the 
minimal UED model studied here, in which the LKP is stable, there have been no direct limits to date. 



2. Experiment 

The Tevatron is a proton-antiproton collider with a center-of-mass energy of ^/s = 1.96 TeV. In Run II, from 
2001 through Sep. 2011, it delivered more than 10 fb~^ to each of its two experiments, DO and CDF. The 
DO detector [10] is a general purpose detector with nearly An coverage. The innermost system is the silicon 
microstrip tracker (SMT), surrounded by the central fiber tracker (CFT) which utilizes scintillating fibers in 
stereo configuration. These systems form the central tracker. Outside of the tracking system is the LAr/Ur 
sampling calorimeter, which is split into central and endcap regions. Beyond the calorimeter is a dedicated 
system for triggering on and detecting muons, consisting of mini drift tubes (MDTs) and proportional drift tubes 
(PDTs). The DO detector utilizes toroidal and solenoidal magnets of each nearly 2 T to curve the trajectories of 
charged particles in the tracker or the muon system, respectively. 



3. Samples 

In this analysis, 7.3 fb^^ of DO data are used, collected in Run II until July 2010. All events which fired a 
single muon trigger and pass data quality requirements are used. 

Most standard model backgrounds are modeled with Monte Carlo. M^+jets, Z+jets and tt events are 
generated with ALPGEN [11] interfaced to pythia [12] for hadronization and showering, while diboson events 
{WW, WZ and ZZ) are simulated using pythia alone. In both cases, the CTEQ6L1 parton distribution 
functions [T^ were used. The W/Z + jets and diboson samples were scaled to a next-to-leading-order (NLO) 
prediction from MCFM [l4l I15j . Cross sections for in next-to-next-to-leading-order (NNLO) were obtained 
from [161. Signal MC was generated for 9 evenly spaced points from i?"^ = 200 GeV to 320 GeV, using pythia 
6.421 and CTEQ5L PDFs [T7]. 

All MC samples are processed through a full GEANT [IB] based detector simulation, and events are reconstructed 
using the same algorithms that are applied to data. Several reweightings are applied to MC to improve agreement 
with data. MC is reweighted so its instantaneous luminosity distribution matches the luminosity profiles of data. 
In the Z-|-jets sample, a correction is applied to the distribution of the Z boson. Finally, reweightings are 
applied to to account for differences in trigger and object identification efficiencies. 

3.1. Preselection 

As a preselection, events with two same-sign muons which fulfill certain reconstruction quality criteria are 
selected. Muons must have a transverse momentum p^ > h GeV/c. Hits in the muon system are matched to 
central tracks, and the global fit of the muon track must have a < 9.5. There must be at least two CFT hits. 
A cosmic veto is applied, where the time between scintillator hits must be less than 10 ns. The distance of 
closest approach (dca) of the muon track to the primary vertex must be less than 0.04 cm, or if there is no SMT 
hit, less than 0.2 cm. Finally, a cut \z{^i) — z{^2)\ < 1 cm is applied, where z is the distance between the muon 
and the primary vertex along the beamline. 

For modeling and suppression of the multijet background, muon isolation is an important quantity. Calorimeter 
isolation IcaX is defined as the sum of Et in the calorimeter in a hollow cone 0.1 < Ai? < 0.4 around the muon, 
divided by the muon p^ -c. Here, Ai? is a measure of angular distance, given by Ai? = \J l\rf + Atp^. Additionally, 
track isolation /tik is defined as the muon p^ of all tracks in the cone Ai? < 0.5 around the muon, divided 
by its p^ . With these isolation variables, the muons are separated into three categories: A tight muon (T) has 
^cai < 0.4 and Jtrk < 0.12. A loose muon (L) is not tight, and has /cai < 0.4 and /trk < 0.25. A non-isolated 
muon (N) is neither tight nor loose. 
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From the preselection, a main selection sample 5" is formed, containing all events where one muon of the 
same-sign pair is tight, and one is tight or loose (TT or TL). Cuts are applied on Acp between the two muons: 
Events with A(/? < 0.25 are removed, and to cut away some of the multijet background, also back-to-back events 
> 2.9). Only events with both muons in the central region Irydotl < 1-5 are used, since the tracking efficiency 
is much lower and less well modeled outside this region. Finally, a minimum cut of > 10 GeV/c is applied. 



4. Background estimation 

While most backgrounds are modeled with MC, for the estimation of the multijet and charge misidentification 
backgrounds, data-driven techniques were used. The modeling of these backgrounds shall be described in the 
following. 

4.1. Multijet 

The QCD multijet background is mostly made up of muons emitted from jets from bb/cc decays. One 
characteristic of this background is that, since jets tend to be balanced, the muon pairs are mostly back to back, 
whereas the angles of signal muons are much less correlated since they come from cascade decays. Since multijet 
muons originate in jets, they tend to be less isolated than signal muons. However, secondary muons from jets 
can become isolated and be misidentified as prompt muons. The probability of this is higher for lower muons, 
and so the multijet background is larger in that region. These properties can be used to model the multijet 
background from data. 

A multijet-enriched sample Q is defined with the same criteria as S, but the following exception: the first 
muon is tight as in the S sample, but the second muon is non- isolated (TN). This sample provides the shape 
of the multijet background, while the normalization is calculated in an orthogonal region where the multijet 
background dominates. This is done by taking the ratio of integrals of the S and Q samples, while flipping the 
Pj. cut for the most isolated muon (lowest /cai): 

N.^j^, . = 0,1,>2 

The resulting normalization factor is strongly dependent on the jet multiplicity i, thus different normalizations 
Ni are determined for zero, one, and two or more jets. These factors are applied as weights to the Q sample (in 
the region > 10 GeV). 

While the assumption that the multijet background is dominant in the Q sample is valid for low p^, there 
is a significant contribution of electroweak (EW) processes, mainly VF-|-jets, in the higher region. Since 
these events are already considered in MC, this leads to double counting. This is solved by estimating this EW 
contamination, and subtracting it from the reweighted Q sample. The EW contamination in Q is estimated by 
applying the Q selection and the Ni weights to EW MC. The additional physical constraint that the multijet 
background goes to zero (and must not be negative) when p^ (^) — cx) is used to derive a normalization factor / 
for the subtracted EW contamination. The final estimation for the background is thus given by: 

BG^ = Q'^'^''^ -N,- f- Q^^ ■ N, + S^^ , i = 0, 1, > 2 . 

Here, the first two terms give the final multijet estimation. Q'^^*-'^ and Q^'~^ denote the Q selection applied 
to data and MC, respectively, and S^'^ is the estimation for non-multijet backgrounds from MC {W/Z+jets, 
diboson, tt). The leading p^ distribution including the multijet background is shown in Fig.[2j 

4.2. Charge mismeasurement 

A sizable background is due to Z ^ M^/-^^ opposite-sign events, where the charge of one muon gets mismeasured. 
This can especially happen at high p^ when the track of the muon is almost straight, or when additional hits 
from other charged particles or noise are present near the track. A small charge mismeasurement probability 
can lead to a large background in the same-sign channel because of the large number of events in the Z peak. 
In principle, this charge flip (CF) background is taken into account by running the Z+jets MC through the 
detector simulation. However, since charge mismeasurement is not well simulated, it is advisable to confirm the 
MC prediction with an independent estimation, which is determined using a data driven method jl9| . 
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Figure 2: Comparison of data and background estimations for the leading muon distribution. 



The idea behind this method is that the DO detector provides two uncorrelated muon charge measurements: 
One "local" measurement from the track in the muon chambers, and one from the matched track in the central 
tracker. Since the accuracy of the central tracker is much larger, this measurement is usually meant when simply 
referred to the "charge" of a reconstructed muon. 

Data events are split into three categories, depending on how often the local and tracker measurements agree: 
AA if the measurements agree for both muons, AD if there is a disagreement for one muon in the pair, and DD 
if both muons have disagreeing local and track charges. The number of muons in each category is then given by: 
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Here, P^^^ = £\oc ■ £ioc is the probability that in a true like-sign event the measurements agree for both muons 
(AA). Since the tracker is assumed to be much more accurate than the muon system, its inefficiency can be 
neglected, and the probability is the efficiency of the local measurement eioc squared. -P^^ is the respective 
probability for an opposite-sign event. Since all considered events are in the like-sign sample, the tracker 
measurement of one of the muons must be wrong, and since it is an AA event, so must be the local measurement. 
It follows -Paa ~ ^loc ■ (1 ~ £ioc)- Similar expressions can be found for the other probabilities. The fraction of 
charge flip events within the like-sign sample is given by 

/flip = ^flip/(A^flip + A^t rue / 5 

and can be obtained by solving the system of equations ([T]). This fraction multiplied by the total number of data 
events yields an estimation of the charge flip background. The obtained number, 161.7 ± 32.4, is in agreement 
with the number of events from the Z — > /i/x background MC, 170.6. From the difference between both, a 
systematic uncertainty for the CF background is calculated. Since that difference is small compared to the 
statistical uncertainty of the CF estimation, both are added in quadrature to arrive at a systematic uncertainty 
of 20.8%. 



5. Signal extraction 

A multivariate analysis (MVA) is used to take advantage of correlations between variables [20]. The following 
variables are uses as inputs to the MVA: 

• Leading and second transverse momentum (p^^ and ) 

• Missing transverse energy {^t) 

• Dimuon invariant mass (Mpair) 
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• Angle between the two muons in the transverse plane (A(/3) 

• Fit of the track matched to the first or second muon (the background contains many events with 
mismeasured tracks) 

• Number of Jets (iVjcts) 

• Product of missing transverse energy and second {^t 'Pti)^ since in the multijet background either one 
of them tends to be low 

• Topological variables Mti and Mt2, calculated as the transverse mass between and one of the muons: 



Mti,2 := \/2^T •PTi,2(l-cosA(p(/T,/ii,2)) 
These provide a good separation between different types of background. 



Some additional cuts are applied before application of the MVA to reduce specific backgrounds or not well- 
modeled regions: J^t > 25 GeV to reduce the multijet background, p^^ = 10 . . . 200 GeV/c, p^^ > 10 GeV/c, a 
cut on the dimuon invariant mass Mpair = 20 . . . 250 GeV/c^, and finally |dcai_2| < 0.05 cm to reduce pile-up 
and badly measured tracks. 
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Figure 3: Some input va riabl es for the multivariate analysis, (a) leading muon p^, (b) second muon p^, (c) Missing 
transverse energy J^t and (d) muon pair invariant mass. The signal is shown for R^^ = 245 GeV. 



A boosted decision tree (BDT) method was chosen as classifier. Individual BDTs were trained for each signal 
point. The output of the BDT for the signal point at = 245 GeV can be seen in Fig. |4] Higher values mean 
more signal-like events. The separation of signal and background is clearly visible, and there is no significant 
excess of data compatible with a signal. 



6. Conclusions 



A search for universal extra dimensions in the dimuon channel using 7.3 fb^^ of DO data has been performed 
in the like-sign dimuon channel. No excess of data over background was observed. A publication of these results 
including limits on the mUED compactification scale is in preparation. 
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Figure 4: Output variable of the BDT for R'^ = 245 GeV. 



Table I: Sources of systematic uncertainties taken into account: 



Systematic 


Applied to 


Uncertainty 


Luminosity 


All 


6.1% 


PDF 


All 


4.0% 


Jet energy scale 


All 


4.1% 


Muon identification 


All 


2.0% 


Tracking efficiency 


All 


1.0% 


Isolation efficiency 


All 


0.5% 


Trigger 


All 


6.0% 


Z/7 cross section 


Z+jets 


3.5% 


Charge flip 


Z+jets 


20.8% 


VK+jets cross section 


W+jets 


8.5% 


Top pair production cross section 


ti 


15.0% 


Diboson cross section 


diboson 


7.0% 


Multijet estimation 


multijet 


40.3% 
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